Engineering many-body quantum dynamics by disorder 
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Going beyond the currently investigated regimes in experiments on quantum transport of ultracold 
atoms in disordered potentials, we predict a crossover between regular and quantum-chaotic dynam- 
ics when varying the strength of disorder. Our spectral approach is based on the Bose-Hubbard 
model describing interacting atoms in deep random potentials. The predicted crossover from local- 
ized to diffusive dynamics depends on the simultaneous presence of interactions and disorder, and 
can be verified in the laboratory by monitoring the evolution of typical experimental initial states. 
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While well-controlled experiments in solid-state sys- 
tems are lacking, the recent advances in atom and quan- 
tum optics allow the experimentalist to study minimal 
models where single-particle dynamics, many-body inter- 
actions, and disorder can be engineered at will. Ultracold 
bosons or fermions loaded into optical lattices, which re- 
alize spatially periodic potentials [l| , are optimally suited 
to study, e.g., quantum transport across disorder poten- 
tials and possible manifestations of Anderson localiza- 
tion in the mean-field regime [3, Moreover, modern 
experiments reach the regime of strong atom-atom corre- 
lations to investigate many-body quantum effects such as 
interaction-driven phase transitions [1, 4] or interaction- 
induced changes of Landau-Zener tunneling rates [1, @]- 

The Bose-Hubbard model well describes ultracold 
bosons in periodic lattices at small fillings (where a mean- 
field theory is obviously bound to fail) and not too shal- 
low lattice depths [H, H 3, Hi ■ A recent study of an open 
Bose-Hubbard system [6 1 furthermore showed that many- 
body interactions lead to similar decay-rate distributions 
as predicted for single-particle transport in disordered 
potentials. More precisely, interactions in a many-body 
system can substitute for disorder in the diffusive regime 
of quantum transport [6] . 

Here we show how to engineer the dynamical prop- 
erties of a many-body Bose-Hubbard system by vary- 
ing the strength of static disorder. We predict that, for 
an intermediate regime of disorder strength, the system 
shows clear signature of global quantum chaos. The lat- 
ter is quantified by spectral measures of q uantum chaos 
i, 0, i, i, m and transport [H [H [iH. Complexity 
arises in our systems from the simultaneous presence of 
atom-atom interactions and disorder. 

We consider a disordered Bose-Hubbard system on a 
ID lattice comprising L sites, defined by the Hamil- 
tonian 
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The operators a!, ai create and destroy bosons at lat- 



tice site £, respectively. The random on-site potentials, 
{e^}, are chosen from a box distribution in [— e/2,e/2]. 
The deterministic Peierls phase 9 in the kinetic term cor- 
responds to a finite (angular) momentum of the lattice 
or, cquivalently, to the presence of a (effective) magnetic 
potential [Tsj . and mathematically to imposing different 
boundary conditions, twisted as opposed to simply peri- 
odic, onto the standard model, Eq. ^ with 9 = p^ . 
Phases could be controlled experimentally as described 
in ^IG'I, while the periodic boundary conditions assumed 
for ID) in the following could be implemented in optical 
ring lattices The effects of {eg} and 9 are in some 

sense complementary. This can be understood switching 
to the Pock basis of the quasimomentum (QM) opera- 
tors 6„ = afe'2'^^«^/-C' diagonalizing the 
kinetic term in ([T|)- In this reciprocal basis, the interac- 
tion term is block-diagonal, the blocks being labeled by 
the total QM= nbl^b^ mod(L) - where the mod oper- 
ation guarantees a QM in the unit interval. Conversely, 
the local potential term, {ee}, couples blocks of different 
QM, since it corresponds to a sum of nonlocal operators 
bl-bfj [l^. Any such operator couples blocks whose total 
QM differs by {k — rf) units, and since b^b,^ is a num- 
ber operator, the coupling within each block reduces to 
a trivial constant. Hence, {e^} induces interblock cou- 
plings, whilst 9 affects the diagonal blocks of fixed QM. 

We are interested in the global dynamics generated by 
([T|)- Our approach to characterize the quantum trans- 
port in the system is twofold: first, we study the spectral 
properties of U) , and secondly, we present results of the 
time-evolution of initial states which are not eigenstates 
of ([T])- The evolution of typical experimental observables, 
such as the spatial population and the mean momentum 
of the condensate particles P, H, 3 , allows one to directly 
probe the here predicted diffusion-localization transition 
(DLT). In contrast to [l^, we do not just focus on the 
regime of small disorder e <C J, and we include 9 to gen- 
eralize the assumed periodic boundary conditions. 

For the case without disorder and filling factors of 
order one, a cross-over between regular and quantum 
chaotic spectra was predicted in [8| when varying the 
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ratio U/J. Quantum chaotic spectra are identified by 
their statistics, more precisely the distribution P{s) of 
the normalized level spacings s = AE/AE follows a 
Wigner-Dyson (WD) distribution Q. Also by adding 
a linear force to the Hamiltonian, a transition between 
regular and chaotic motion can be identified by a cross- 
over from Poisson to WD statistics 0, 0] , corresponding 
to a regime of strong (Stark) localization or of quantum 
chaos, respectively. Based on translational invariance, 
most previous results [1, 0, [l] concern a subset of states 
with constant (conserved) QM. Since in the experiment 
it is generally hard to focus on a subset of states tak- 
ing part in the dynamical evolution, it is desirable and 
more general to search for global quantum chaos, which is 
not just restricted to independent subsets of the system's 
eigenvalue spectrum, e.g. corresponding to constant QM. 
Disorder (which cannot be completely avoided in any real 
system) naturally breaks the translation invariance [l^ . 
and we will show that the dynamics of a dilute boson 
system (induced by its spectral properties) can be con- 
trolled by the combined action of atom-atom interactions 
and the random potential in P]). 

For a complex spreading of the system's eigenf unctions 
in the eigenbases of the integrable cases ( J = = e) and 
{U = = e), the energy scales defined by the terms in 
P]) should be roughly of the same order of magnitude, 
i.e., J ~ f7 ~ e. For such a situation and small filling, 
we indeed found clear signatures of global quantum chaos 
in the system. Our results are shown in Fig. [U which 
collects avoided crossing scenarios of the energy levels (a) , 
the cumulative distribution C(s) = ds' P{s'), and the 
number variance of the energy levels (which measures 
the long-range correlations in the spectrum _10] ) in (b,c). 
Both C{s) and E^ well agree with the WD predictions 
for a Gaussian Orthogonal Ensemble of random matrices 
[?,, 8, .9, .10]. Residual symmetries of ([T]) (e.g., a reflection 
symmetry for special values of the number of atoms N 
and the system size i fs']) can be destroyed hy 9 ^ 0. 

Our results of Figs. [T] and O depend on the procedure 
chosen to unfold the energy spectrum (i.e., to make the 
density of states approximately constant), which is neces- 
saryto compare to the normalized theoretical predictions 
d, [l3| ■ We used a rescaling of the levels by the numer- 
ically obtained local density of the raw data, which is 
independent of further assumptions on the original level 
density. Choosing a small window (over 5 ... 10 levels) for 
computing the local average permits an optimal match on 
smaller scales, whilst larger windows (over 20 ... 40 lev- 
els) are chosen to compute the number variance. Our 
analysis considered all levels of a given spectrum at fixed 
parameters, and we checked that excluding levels at the 
band edges does not change our results, as long as we 
stay in the dilute limit of UN/L < 1. 

For e < 0.5 and e ^ 1, we observed a trend towards 
globally regular dynamics, a consequence of good, yet not 
perfect Poisson statistics. In the limit e 0, the various 
QM blocks uncouple, and as exercised in @, 0, Q one has 
to concentrate just on one of these blocks in the spectral 
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FIG. 1: (color online) (a) Sector of the spectrum in the band 
center as a function of U, for A'^ = 3, L = 15, J = 1, ~ 0.119 
for e = 1 (a linear function of U was added to E to eliminate 
an overall trend), (b) C(s) (obtained from collecting levels 
of 25 disorder realizations) and (for one realization in the 
inset) at (/ = 1, corresponding to the WD prediction (solid 
lines), (c) same as in b) but for a regular case with TV — 
3,L = 15 (o), 20 (o), (7 = 1 = J and e = 10,61 = 0, together 
with the Poisson predictions (solid lines). 



analysis. Any small e 7^ 0, however, destroys the transla- 
tion invariance, making an analysis of the full spectrum 
quite intricate. We, therefore, concentrate for a moment 
on the case of large disorder. As shown in Fig. [TJc), this 
limit is well characterized by a Poisson distribution. This 
result is expected, since the eigenstates become pinned at 
the randomly distributed minima of the potential, lead- 
ing to a small residual overlap between them. Of course, 
for finite U and L, such a localization cannot be per- 
fect [l3|- Indeed, we observe a better correspondence in 
our system with the Poisson prediction (see inset of Fig. 
[ijc)) for smaller fillings, consistent with single-particle 
localization theory J(]J. This trend was observed when 
increasing L < 20 at fixed N — 3 and 4, or when decreas- 
ing 1 < iV < 4 at fixed 10 < L < 20. 

The DLT occurs for all data sets shown in Fig. [5] at a 
critical value ~ 4. We characterize this crossover by 
a statistical test [l^, which measures the deviation 
from WD. The saturation of our Xwu measure for large 
e, at the right of the graphs in Fig. ^ signals the con- 
vergence to a Poisson distribution, an example of which 
is explicitly shown in Fig. [BJc). For fixed N and L, ecr 
scales linearly with U = J, as expected since our Hamil- 
tonian in ([T]) is scale invariant for a given realization of 
disorder (i.e., U = J defines the energy scale for e). At 
fixed U = J, the crossover to the localized regime de- 
pends on the filling, and our results from Fig. [2] suggest 
the following functional form XwDi^)^ with x = (.\J L/N^ 
in the range e ~ 4. . . 15. It is numerically hard to ob- 
tain a full scaling function for interacting systems with 
larger N and L, since it is necessary to diagonalize the 
full system, not just one QM block as in [1, 0, Q. 




0.8 

□ o 

-§ < 0.2 

i 


ft 


— ■ — Berry-Robnik 
^ Semi-Poisson 

12 3 4 

S o 


4 ^ 

o <S 
o, Bo-° 

' ° ^"■^ 
JO, 

(b) 



10 20 30 40 50 

FIG. 2; (color online) (a,b) statistical test (with val- 
ues close to zero for good WD statistics) as a function of the 
scaled disorder parameter, for (a) L = 15 and A'^ = 2 (□), 3 
(o), 4 (A), and (b) for A*- = 4 and L = 8 (□), 10 (o), 14 (a), 
17 {<>), at J = 1 = U,9 = 0. The inset in (a) shows a zoom 
of iV = 3 for small €,9 = (o) and 61 ~ 0.119 (+). Each data 
point is averaged over 20 disorder realizations. The inset in 
(b) shows the relative deviation for the marked point in the 
crossover regime (dashed line) from two heuristic interpolat- 
ing laws between WD and Poisson: Semi-Poisson (grey) [2ll ] 
and Berry-Robnik 22] (black). As expected [2ll. [23] ■ corre- 
spondence is not perfect with neither of those laws, but the 
Berry-Robnik lies closer overall. 

Interestingly, our spectral analysis of Figs.[T]and[2]does 
not show a dependence on 6 in the localized regime (e > 
4), whilst especially for e < 1 (see inset of Fig.[2][a)) both 
couplings by ^ and by e 7^ can conspire to enhance 
the quantum chaotic properties of the full spectrum (i.e., 
not only of a subblock of fixed QM) . We therefore can use, 
to some extent, both parameters as independent handles 
to change the global spectral properties. 

The dependence of the spectrum on the choice of 
the boundary conditions defined by 9 is reminiscent of 
the Thouless conductance, another prominent measure 
to characterize the transition between extended and lo- 
calized states [ill. We computed the Thouless con- 
ductance, which essentially is given by the curvature 
Ct = {\(fE/de^\) « {\2[E{e) ~ S(0)] 1/02^^0^ 9 ^ Q 
prl [ist , which was geometrically averaged [ll| over the 
full spectrum and 40 realizations of disorder. Our results 
are shown in Fig. [3l In contrast to the distributions of 
nearest level spacings, where the small but finite mixing 
of QM blocks at e ^ does not allow to well character- 
ize the true type of dynamical regime, the curvature is 
a local property of the spectrum. Hence, in the diffusive 
limit e ^ 0, we find the expected d 
with a « 2...1.8 (A^ = 2...4,L = 15) and a « 1.3 
(A^ — A:,L — 10). In analogy to Fig. [21 the crossover 
between the diffusive (quantum chaotic) and the local- 
ized regime sets in at Ccr ~ 4. For e > 4, our results 
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FIG. 3: (color online) Level curvature Ct as a function of 
the disorder strength for A^ = 2 (□), 3 (o), 4 (A) and 
J = 1 = U,8 = 0, L = 15. Two different scalings are cho- 
sen to highlight the DLT around the dashed line. Increasing 
the filling (particularly for N = 4:,L = 10, shown as full di- 
amonds) leads to deviations from the single-particle scalings 
for both small and large disorder. 

confirm an exponential scaling (typical of finite-size lo- 
calized systems) Ct oc exp(— const • e^), with (3 « 0.8 
(A^ = 2 . . . 4, L = 15) and /3 « 0.5 (A^ 4, L = 10). The 
systematic deviation of both exponents from the single- 
particle predictions {a — 2 and (3=1 [lH) with increas- 
ing filling factor highlights the structural change of the 
level dynamics in the presence of atom-atom interactions. 
In the following, we focus on experimentally observable 
consequences of the spectral analysis presented so far. 
Fig. [4| shows the temporal evolution of initial states, typ- 
ically prepared in experiments with Bose-Einstein con- 
densates. Panel (a) presents the mean momentum on 
the lattice 0, [H, , defined as 

Pit) = ^m)\ E - ^-c-) i^w) 

in the direct and reciprocal space, respectively, for 
\i; {t = )) = (6t)^lO)/\/M, with n = 2. For strong 
disorder, the momentum decays almost instantaneously 
to zero, and is further characterized by small, random 
fluctuations. The quantum chaotic behavior for e = 1 is 
visible in the strongly correlated large-scale fluctuations, 
characterized by a slowly, algebraically decaying Fourier 
transform of the time series p{t), c.f. Fig. [H^b). The lat- 
ter implies a large number of modes being present in the 
evolution oip{t), a standard signature of complex dynam- 
ics ^2^. Complex transport behavior was predicted also 
in 25] by analyzing the power spectrum of oscillations 
in a three-well system. In contrast to (a,b). Figs. [Hc,d) 
show the real space dynamics of a box-distributed initial 
state with one atom in wells £ = 7, 8, 9 and none else- 
where. Whilst the chaotic mixing of all wells dominates 
for e = 2, in the localized regime e = 10 we observe the 
expected strong pinning at the random minima of the 
potential (particularly in the 7th well in Fig.[3I^d)). The 
fluctuations in the latter case arise from our finite values 
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FIG. 4: (color online) (a) p{t) for a three particle initial 
state with fixed QM, for L = 15, J = 1 = ?7, 6^ = 0, and (b) 
Fourier transforms of the curves from (a), for e = 1 (thick 
lines) and e = 10 (thin grey lines). The evolution of the well 
populations for a box initial state of = 3 atoms in real 
space is shown in (c) for e = 2 and in (d) for e = 10 (same 
parameters as in (a), color code defined in (e)). The time unit 
is r = 2-Kh/ J[Eb\ ~ 15 ms for a lattice constant of the 
order 412 nm and a lattice depth of 10 recoil energies Er for 
Rb atoms P, I3- The inset in (a) highlights the slow decay 
for e = 1. 



turn chaos (without any residual symmetry in the system) 
and strong localization. Yet, the crossover between the 
two regimes is systematic, and as exemplified, even the 
evolution of single, reduced experimental observables of 
our many-body problem can be used to directly visualize 
the change in the quantum spectra analyzed in Figs.[T][31 
Of experimental relevance is in particular the difference 
in the short-time evolution of p{t) for experimental de- 
tection times < 1 s. 

In summary, we showed how one can scan between the 
different dynamical regimes of the Bose-Hubbard system, 
characterized by global quantum chaos and by essentially 
localized bosons, by varying the strength of static disor- 
der. The dynamics of initial states which are far from 
eigenstates of the system could be used as a clear ex- 
perimental signature of this crossover. Moreover, our 
analysis of the many-body level curvatures opens a new 
link to transport problems in mesoscopic solids [111, [l^ 
and photonic lattices [26^ , where global chaotic properties 
are accessible by conductance measurements. As adum- 
brated in the here presented, extended and unify- 
ing characterization of the spectral properties of a disor- 
dered many-body problem, may be useful to obtain, for 
instance, experimentally accessible estimates for the lo- 
calization properties for such com plex systems. However, 
the necessary scaling arguments flzl as a function of N 
and L make such an approach challenging for up-to-date 
computational resources. 



of U and L. For a typical realization of {e^}, the evolu- 
tion will be asymmetric as seen in (c), which allows one 
to distinguish it from the e = case. Equivalently, whilst 
for e = QM is conserved, for e ^ QM starts to devi- 
ate from its initial value, which is for the data in (c,d). 
Fig. m presents the limits of fully developed global quan- 
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